December 23, 2009 15:2 World Scientific Review Volume - 9in x 6in Spin-Orbit-Book 



ON 

o 
o 

(N 

o 

Q 

m ! Chapter 1 

(N 

Spin-Hall Effect in Chiral Electron Systems: from 
^ ■ Semiconductor Heterostructures to Topological Insulators 

Q , P.G.Silvestrov 1 and E.G.Mishchcnko 2 

1 Theoretische Physik III, Ruhr-Universitat Bochum, 44780 Bochum, 

Germany 

2 Department of Physics, University of Utah, Salt Lake City, Utah 84112, 

USA 

<— I 1 The phenomenon of mesoscopic Spin-Hall effect reveals in a nonequilib- 

Q ' rium spin accumulation (driven by electric current) at the edges of a 

O ■ ballistic conductor or, more generally, in the regions with varying elec- 

tron density. In this paper we review our recent results on spin ac- 
cumulation in ballistic two-dimensional semiconductor heterostructures 
. with Rashba/Dresselhaus spin orbit interactions, and extend the method 

QQ ' developed previously to predict the existince of spin-Hall effect on the 

, surface of three-dimensional topological insulators. The major differ- 

■ ence of the new Spin-Hall effect is its magnitude, which is predicted to 

be much stronger than in semiconductor heterostructures. This happens 
^v^j , because in semiconductors the spin accumulation appears due to a small 

spin-orbit interaction, while the spin-orbit constitutes a leading term in 
Q\ ' the Hamiltonian of topological insulator. 

o ' 

' 1.1. Chiral electron systems 

b . 

Chiral systems feature electron states whose quantum properties (e.g. spin) 
depend on the direction of propagation. One example of such chiral states 
arises from spin-orbit coupling that originates from relativistic corrections 
to the dispersion law of band electrons. In particular, for a typical two- 
dimensional electron gas (2DEG) the intrinsic asymmetry of a confining 
quantum well geometry is accompanied by a strong perpendicular (z) elec- 
tric field that leads to spin-orbit coupling of the Rashba type, 1 

H R = AZ • (<T X p) = X(a x Py - (JyPx). (1.1) 
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As a result, the spin degeneracy is lifted via effective momentum- dependent 
Zeeman field, h p = \(—p y ,p x ). The electron spin eigenstate is thus deter- 
mined by its momentum, p. 

Another important type of spin-orbit coupling occurs in 2DEG formed 
by semiconductors with broken inversion symmetry, e.g. GaAs, InAs. While 
only third order in momentum in the bulk GaAs, this coupling, known as 
Dresselhaus interaction, 2 is "upgraded" in two dimensions to the linear 
order by virtue of transverse momentum quantization and becomes, 

H D = ^d(o~xPx - O-yPy)- (1-2) 

Chiralities acquired from different spin-orbit couplings and (|1.2|) are 
opposite, in a sense that the electron wave function acquires opposite Berry 
phases, ±7r, upon (adiabatic) completion of a loop in the momentum space 
(enclosing the degeneracy point p = 0), depending on whether a > (3 or 
vice versa. 

Dirac fermions in graphene represent another realization of the chiral 
system. 3 Due to hexagonal symmetry of underlying two-dimensional hon- 
eycomb atom arrangement the low-energy electron excitations are combined 
into two Dirac cones (valleys), K and K, within the first Brillouin zone. 
The corresponding effective Hamiltonians are 

H K = v(r x p x + Typ y ), = v(r x p x - T y p y ). (1.3) 

Here r stands for the pseudospin operator that acts in the sublattice space. 

It is interesting to note the analogy between total graphene Hamiltonian 
and a semiconductor with equal Rashba and Dresselhaus couplings (a = (3) 
in the case of a conventional 2DEG. The latter case features spin eigenstates 
that are momentum-independent 4 and, thus, somewhat trivial. However, 
since the two cones in graphene are separated in momentum space, the chiral 
physics can still be observed. (Yet some phenomena are trivially absent, 
e.g. pseudospin-Hall effect, cf. last section). It is important to emphasize 
that the chirality in graphene has nothing to do with spin-orbital coupling 
(which is rather weak in carbon allotropes) and is a consequence of the 
crystal geometry. 

Another notable example of a chiral electron system is a topological 
insulator. The latter is different from the usual band insulator in that its 
valence and conduction bands originate not from different atomic orbitals 
but from the same spin- orbit- split orbital. In 2D HgTe quantum wells this 
gives rise to topologically protected edge states leading to a recently pre- 
dicted 5 and discovered 6 quantum spin-Hall effect. Yet even more intriguing 
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twist has recently been added by a discovery of two-dimensional states on 
the surface of 3D topological insulators 7 ' 8 Bio.gSbo.i, Bi2Sb3, and Bi2Te3. 
The fundamental difference from graphene is that the number of pockets of 
the Fermi surface within the Brillouin zone is odd, with the complementary 
species residing on the opposite surface of a sample. Thus, these states are 
not simply separated in the momentum space (like in graphene), but also 
separated in real space. This removes the above mentioned "trivialization" 
that is present in graphene. In particular, spin-Hall effect can occur in 
topological insulators. As confirmed by the first-principle band structure 
calculations 9 spin structure of these states is indeed chiral. In particu- 
lar, the low-energy Hamiltonian of Bi2Te3 can be deduced from symmetry 
considerations to be of the form, 46 



For an ungated and undoped Bi2Te3 the last term is generally of the same 
order of magnitude as the first two. Still, as a starting point it is useful to 
neglect the effects of anisotropy. 

1.2. Spin-Hall effect 

Spin-Hall effect 10 is the name given to a class of phenomena that ex- 
hibit boundary (surface, edge) spin polarization when electric current flows 
through a system with significant spin-orbit interaction. It has been ob- 
served in both 3D 11-13 and 2D systems. 14 It is customary to distinguish 
two mechanisms that could lead to this effect. The extrinsic mechanism 
is the dominant one in 3D semiconductors and originates from scattering 
off impurities. 15-18 Presence of impurities is unavoidable in high carrier 
density 3D semiconductors that rely on doping. Quite contrary, intrinsic 
mechanism 19,20 originates from spin-orbit-split band-structure, which in- 
duces spin precession when electric current is driven through the system. 
This mechanism can in principle be realized in coexistence with ballistic 
transport in 2D electron systems systems. Indeed, by placing dopant far 
enough from GaAs/AlGaAs interface one can reduce effects of disorder 
scattering. 

In the present paper we concentrate on the intrinsic spin-Hall mecha- 
nism. It is known, however, in two-dimensional electron systems with spin- 
orbit coupling linear in momentum (typical for n-doped heterostructures) 
any scattering that leads to a stationary electric current via deceleration 
of electrons by impurities, phonons, etc., will negate the precession due to 



H T i = v(o- x p y - <j y p x ) + 




2m 



+ a<T z p x (p 2 x + 3p 2 y ). 



(1.4) 
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external electric field and result in the exact cancellation 21-25 of the bulk 
spin-current in a dc casj^. 

There are several ways to avoid such cancellation, in particular, to use 
ac currents with frequencies exceeding the inverse spin relaxation time. 21,26 
The second possibility and the one of interest to us here, is make a system 
sufficiently small and clean (ballistic) so that the electron mean free time 
exceeds the time of flight across the systems. The corresponding scenario is 
known as the mesoscopic spin-Hall effect. 27 While initial theories of spin- 
Hall effect in infinite systems had addressed the auxiliary quantity of spin 
current (for a review see Refs. 28,29 ), in a finite geometry it is both easier 
and more relevant to calculate spin polarization density, which is an exper- 
imentally measurable quantity. 30 Such edge polarization was considered 
by numerical methods in several earlier publications 27,31-34 as well as both 
analytically and numerically in our previous papers. 35,36 

A crucial note is due. The edge spin polarization in ballistic systems 
appears not as a result of electric field-driven acceleration of electrons and 
associated with it precession of spins. Indeed, electric field in a ballistic 
system is absent as the electric potential drop occurs over a contact region 
with the leads rather than over a bulk of a ballistic conductor. Nevertheless, 
spin precession does occur. It originates from accelerated electron motion 
in the boundary potential that provides lateral confinement. Bias applied 
to the contacts ensures that populations of left- and right-moving states are 
different and the net spin precession appears. Naturally, it is proportional 
to the applied bias V. The net spin accumulation near an edge of a 2DEG 
is independent of the shape of the boundary potential, 



where vf is the bulk value of the Fermi velocitjo Spin accumulation ap- 
pears in the second order in spin-orbit interaction and vanishes for equal 
Rashba and Dresselhaus coupling strengths. 

Position-resolved spin polarization s z {x) can be found analytically in 
two important situations. First is the case of a smooth confining potential 
U(x) and the second is an infinite hard- wall boundary. We begin with 

a This exact cancellation does not occur in 3D- or 2D hole-systems that feature non-linear 
spin-orbit couplings. 

b Strictly speaking, the integration in Eq. Ill.5|l goes from a pint far outside the 2DEG 
(— oo) to some point deep inside 2DEG (+00), but still far away from the edges. Inte- 
gration across the whole conductor would give zero, reflecting the fact that spin accu- 
mulation at the opposite edges has opposite signs. 




(1.5) 
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analyzing semiclassical electron motion in smooth potentials in Sec. 11.31 
In Sec. 11.41 we derive the result (|1.5|) for the net spin accumulation and 
illustrate it using the case of hard-wall boundary. In Sec. 11.51 we present 
a method of kinetic equation that allows to find local spin-polarization for 
smooth boundary potentials. In Sec. 11.61 singular dynamics near classical 
turning points is discussed. Finally, in Sec. ll.7l we use methods developed in 
the preceding chapters to describe spin-Hall effect in topological insulators. 

Many of the results presented in this paper were published in the journal 
articles Refs. [35-37]. However, the derivation of the spin accumulation via 
the Kinetic equation in Sec. 11.51 and prediction of the nonequilibrium spin 
accumulation in topological insulators in Sec. 11.71 are presented here for the 
first time 

1.3. Semiclassical electron motion 

Consider gated 2DEG with Rashba spin-orbit interactional, described by 
the Hamiltonian 

H = 1- \(p y a x -p x (J y ) + — ^— +U(x,y). (1.6) 

Potential U(x, y) is created by the external gates, or the edge potential 
ensuring the in-plane confinement of 2DEG. Possible effect of disorder on 
U(x, y) are going to be neglected. 

Classical electron dynamics described by the short wave length limit of 
the Hamiltonian (|1.6|) reveals a number of very unusual features. As a first 
step, in this section we demonstrate how the trajectories-based approach 
allows to describe propagation of fully (in-plane) polarized electric currents 
through mesoscopic constrictions. 

Construction of semiclassical solutions of the Schrodinger equation 
with the Hamiltonian (jl.6[) follows the reasoning of the conventional 
WKB approach, 37-41 which is valid for a smooth potential, ?i|VE/| <§; 
min(p 3 /m,p 2 A). 

Without the external potential U, the electron spectrum consists of 
the two subbands, E±(p x ,p y ) — (p ± mX) 2 /2m. The subbands meet at 
only one point, p = 0, and the spin in each subband is always aligned 
with one of the in-plane directions perpendicular to the momentum p. The 
semiclassical electron dynamics 38 naturally captures the essential features 
of this translationally invariant limit. The classical motion in each subband 

c Similar methods may be used to describe more complicated interactions, like Dressel- 
haus, 2 or combined Rashba and Dresselhaus interactions. 
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is determined by the equations of motion which follow from the effective 
Hamiltonian: 

H cS = ^l + U( X ,y). (1.7) 
Zm 

Despite the fact that spin does not appear in this equation, one can easily 
construct semiclassical wave functions, which have spin pointed within the 
xy plane perpendicular to the momentum: 

iP = ue iS /\ u = M( y% +5L V (1.8) 
V 2 P \ ± VPy ~ l P* J 

Here the action S is related to the momentum by p — VS 1 , and p = u^u 
is the classical density for a family of classical trajectories corresponding 
to a given energy E. The action S obeys the classical Hamilton- Jacobi 
equation, 

(|Vff|±mA) 2 

+ U{x,y)=E. (1.9 

zro 

During its motion, an electron described by Eq. (|1.8p changes the mo- 
mentum p but always remains in the same spin-subband. To change the 
subband the electron trajectory should pass through the degeneracy point 
where both components of momentum vanish simultaneously, p = 0, which 
is generically impossible. Moreover, with the proper use of potential bar- 
riers, one may realize a situation where electrons of only one subband are 
transmitted and the others are totally reflected. This leads to strong in- 
plane polarization of the transmitted electron flow. 

To take into account the out-of-plain spin precession one has to go 
beyond the approximation Eq. (|1.8p . as was done in Ref. [37]. Instead 
of doing so we will describe in Sec. 11.41 a method allowing to find easily 
the expectation value of a z for potential depending only on one coordinate 
U = U(x) (boundary potential). 

1.3.1. Sharvin conductance 

To give an example of a spin-polarized current let us consider transmission 
through a barrier, U{x), varying along the direction of a current prop- 
agation. We assume periodic boundary conditions in the perpendicular 
direction (y + L = y), thus p y is an integral of motion^. For a smooth 

d Spin polarized currents on a cylinder with x-dependent spin-orbit interaction were 
considered recently in Ref. [42] 
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Fig. 1.1. Conductance (in units of e 2 /h), and spin polarization of the current vs gate 
voltage (in units of mA 2 /2). Dashed lines show the smoothed curves 111.141 1. (11.151 1. solid 
lines show the quantized values for mXL/h = 10.5-7T. Dotted line shows the conductance 
without spin-orbit interaction. 

potential U(x) the conduction channels may cither be perfectly transmit- 
ting or completely closed. The conserved transverse momentum takes the 
quantized values, = 2Tthn/L. Consider the functions 

El(p x )= (pn *™ X)2 , P n = ^/pl+pf (1.10) 

For « / the function E±(p x ) splits into two distinct branches. At any 
point x the equation 

El( Px )=E F ~U(x) (1.11) 

yields solutions p x and p x , corresponding to left- and right-moving elec- 
trons. Application of a small bias implies, e.g., the excess of right movers 
over left movers far to the left from the barrier. Particles are transmit- 
ted freely above the barrier if Eq. (|1.11| has a solution, p x , for any x. 
Let fi = Ep — JTmax be the difference between the Fermi energy and the 
maximum of the potential. The nth channel in the upper branch opens 
when 

fi = (2tt%| +m\L) 2 /2mL 2 . (1.12) 

For the lower branch E"(p x ) Eq. (|1.11[) has four solutions (two for right and 
two for left movers) for \n\ < mAL/2irh and x close to the top of the barrier. 
However, far from the barrier (where the excess of right-movers is created) 
there are still only two crossings described by Eq. (|1.11[) , one for right and 
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one for left movers. As a result, all the extra electrons injected at x = — oo 
follow the evolution of a solution of Eq. with the largest positive p x . 

For all |n| < mXL jlixti such a solution does exist for any positive /x. Thus, 
at fi = as many as no = mXL/irh channels open up simultaneously. The 
channels with higher values \n\ > m\L/2irh in the lower subband open 
when 

H = (2irh\n\ - m\L) 2 /2mL 2 . (1.13) 

According to the Landauer formula, ballistic conductance is given by the 
total number of open channels multiplied by the conductance quantum 
G = e 2 /h 

r -r —i V^rn + mX, < p < mX 2 /2 

U ~ U %h\ 2V2^, p>mX 2 /2. U ' 14J 

This dependence G{p) is shown in Fig. 11.11 The striking evidence of the 
presence of spin-orbit interaction is the huge jump of the conductance at 
the pinch-off point, as opposed to the conventional square-root increase 
in the absence of spin-orbit coupling. This jump is a consequence of the 
"Mexican-hat" shape of the spectrum E_(p x ,p y ). Accuracy of Eqs. (|1.12|) 
and (|1.13p is sufficient to resolve the steps in the conductance due to the 
discrete values of |n| = 0, 1,2, (conductance quantization), as shown in 
Fig. 11.11 The steps in G(/i) are abrupt in the limit dU jdx — > 0. 

Close to the pinch-off, at /i < mX 2 , the conserved p y component of 
the electronic momentum varies for different transmitted channels within 
the range \p y \ < mX. Therefore, far from the barrier, where the Fermi 
momentum is large pf S> mX, we have p x 3> p y and transmitted electrons 
propagate in a very narrow angle interval \9\ < \J mX 2 /2Ef <§C 1. Since 
the electron spin is perpendicular to its momentum, we conclude that the 
current due to electrons from each of the subbands is almost fully polarized. 
The total polarization of the transmitted current is given by the difference 
of two currents 

= W<r y v x 1>) = min(l! y^V2^), (1.15) 

{tp^v x ip} 

which is also depicted in Fig. II. II This current polarization may also be 
viewed as a creation of in-plain noncquilibrium spin density, maximal on 
the barrier. 

Derivation of Eqs. (| 1 . 14[) and (| 1 . 1 5[) was greatly simplified because of the 
periodic boundaries. Our next example demonstrates semiclassical treat- 
ment of realistic boundary conditions. 
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1.3.2. Quantum Point Contact 

Let us consider probably the most experimentally relevant example of a 
quantum point contact, described by the potential 

mfl 2 x 2 muj 2 y 2 
U(x,y) = — + — - — . (1.16) 

Even in this simple model the electron flow in the presence of spin-orbit 
interaction acquires a number of interesting and peculiar features. This 
become clear already from the figure \1.2\ where we show three types of 
trajectories in such potential. Each kind of trajectories is necessary for 
calculation of conductance. 

Classical equations of motion follow in the usual manner from the ef- 
fective Hamiltonian (|1.7|) : f = dH^/dp, p = —dH c s/dr. We consider 
quantum point contact (QPC) close to the opening with only the lower 
_E_ subband contributing to the conductance. A crucial property of the 
Hamiltonian H c g , Eq. (II. 7|) . is the existence of a circle of minima of the 
kinetic energy at \p\ — mX. Expanding around a point on this circle, 
p Xa = mAcosa, p yo = mAsina, one readily finds the equations of motion 
for V = p x cos a + p y sin a — mX <C mX, 

V + (-£l 2 cosa 2 + u 2 s\na 2 )T = , a = 0. (1.17) 

The trajectory is found from the relations, x = V cos a/m , y = V sin aim. 
We observe from Eq. (|1.17[) that only the trajectories within the angle 

tan \a\ < tan ao — il/uj (1-18) 

are transmitted through QPC. Trajectories with larger angles are trapped 
(oscillate) within the point contact. Examples of both types of trajectories 
are presented in Fig. 11.21 Quantization of trapped trajectories would give 
rise to a set of (extremely) narrow resonances in the conductance, specific 
for spin-orbit interaction. Below we consider only the smoothed conduc- 
tance. 

To calculate the current J through QPC one has to integrate over the 
phase space of the states which are transmitted from left to right, 

J =J dy J ev *j£w= GV > ( L19 > 

and have the energy within the interval fi — eV/2 < E- < fi + eV/2, with V 
standing for the applied voltage. In this section we define /i as the difference 
between the Fermi energy and the value of the potential at the saddle point 
/i = Ep — J7(0, 0). The integral is most simply evaluated at x — (with 
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X 



Fig. 1.2. Three kinds of trajectories in the point contact, a, transmitted trajectory 
whose momentum is always collinear with the velocity, b, trajectory bouncing inside 
the QPC. This trajectory is periodic in the linearized approximation described in the 
text, while the exact calculation for finite amplitude shows its slow drift, c, transmitted 
trajectory whose momentum inside the contact is opposite to the velocity. Electrons flow 
from left to right. Arrows show momentum and spin orientations. Few cquipotential 
lines are also shown. 

the velocity given by v x = V cos a/m). The allowed absolute values of the 
momentum are 

2/i - eV - mu?y 2 < V 2 /m < 2^ + eV - mu 2 y 2 . (1.20) 

The angle interval of transmitting trajectories consists of two domains: 
\a\ < a , V > 0, and \a — n\ < a , V < 0. The appearance of the latter 
range of integration is highly non-trivial. A simple reasoning shows that the 
particles with the velocity antiparallel to the momentum (v x > 0, p x < 0) 
should not contribute to the conduction in the case of a transition through 
a one-dimensional barrier U = U(x), see Eq. (|1.14|) . Despite corresponding 
to the right-moving electrons, these states do not originate in the left lead. 
Indeed, they exist only in the vicinity of x = 0, but disappear as x — > — oo 
and, thus, cannot be populated by the excess electrons (except due to 
the tunneling transitions which arc irrelevant in the scmiclassical regime). 
Such trajectories, however, do exist in QPC, Eq. (|1.16|) . as demonstrated 
in Fig. 11.21 After passing through QPC the trajectory bounces at the 
wall reversing its velocity. This kind of classical turning points, where 
both components of the velocity vanish simultaneously, are specific for the 
effective Hamiltonian (|1.7[) . The existence of transmitting trajectories with 
|ct — 7r| < ao, g < results in the doubling of the conductance. Simple 
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calculation yields 

G = G 4mX f nao M. (1.21) 
Trnuj V m 

The presence of a threshold angle an,, as well as the square- root dependence 
of G(fi), are in a sharp contrast to the well-known result G = Go[i/irhu>, i n 
the absence of spin-orbit interaction. 

Since Eq. (|1.17p describes only the linearized electron dynamics, 
Eq. (jl.2ip is formally valid if /i <C mA 2 @ Nevertheless, the current re- 
mains totally polarized for < \i < mA 2 /2 [similar to Eq. (jl.151) ] 

(Vy) = -7TJ 77- = 1- K 1 - 22 ) 

With increasing the chemical potential, /i > toA 2 /2, transmission via the 
upper subband E+ kicks in and the degree of polarization gradually de- 
creases, like it happened in Eq. (|1.15|) . 

In InAs-based heterostructures, typical value of spin-orbit coupling 43 is 
\h = 2x 10 _11 eVYn. Characteristic spin-orbit length In = h/m*\ — 100 nm 
and energy m*A 2 /2 = 0.1 meV. In order to have strongly spin-polarizing 
QPC, the latter should support many transmitting channels at chemical 
potential \x ~ m*A 2 /2 3> fko. This condition can, equivalently, be written 
in terms of the width of the point contact Ay, see Eq. (|1.16)1 . as Ay 3> Ir- 



1.4. Edge spin accumulation. Hard wall 

Consider a semi-infinite ballistic 2DEG described by the Hamiltonian 
Eq. (|1.6p . where potential depending only on one coordinate, U = U(x), en- 
sures boundary confinement (see Fig. II. 3|) . The system is attached to two 
ideal reflectionless leads injecting equilibrium electrons into 2DEG. The 
chemical potentials of the leads are shifted by the applied voltage, eV. The 
current flow along the y-direction in case of spin-orbit interaction Eq. (|1.6p 
results in an edge spin accumulation s z (x). In this section we first derive 
an exact formula for the total amount of spin accumulated at the edge, 
J s z dx. Then we consider in more details the case of hard wall potential, 
U(x < 0) = co, U(x > 0) = 0, where exact results for s z (x) are available^. 



e Still the number of open channels should be large for semiclassics. 

'Nonequilibrium edge spin accumulation at the hard wall boundary of disordered con- 
ductor was considered in Ref. [44] 
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(A) 

*U(x 



(B) 




2DEG | 



E F +eV 




Fig. 1.3. Profile of a boundary potential U(x). B) Geometry of the system: two- 
dimensional electron gas (x > 0) is attached to two ideal reflectionless metallic leads 
filled by equilibrium electrons up to different chemical potentials. C) Spin-orbit-split 
subbands structure. D) Difference in population of left- and right-moving electrons due 
to the applied bias eV. 

1.4.1. Net spin accumulation 

For the potential independent on y (U — U{x)) the corresponding mo- 
mentum component, k y , is an integral of motion. It is convenient to 
use the Fourier representation along the y-axis for the electron operators, 
tfj(r) = J2k Ck y (x)e tkyV . We employ here the second quantization formal- 
ism. One can derive the equation of motion for the expectation value of the 
electron spin operator, s(k y ,x) — \{c\ (x)ack y (x)}, which can be readily 
written in the form, 



Here JJf stands for the conventional operator of spin-current, i.e., 



In a steady state the lhs of Eq. (|1.23|) vanishes. Integrating Eq. (|1.23|) over 
the ^-direction, we obtain for the net spin polarization, 



(1.23) 




"y 



(1.24) 
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It is straightforward to calculate the value of the (fc y -resolved) spin current 
J%(k y , oo) inside the bulk of a 2D system: 




where np(k x ,k y ) stands for the population of different momentum states 
in the subband /3. Only "uncompensated" states contribute to the non- 
equilibrium spin polarization given by Eqs. (|1. 24111.25)) ; these states describe 
electrons that originate in the left lead (k y > 0) and belong to the energy 
interval near the Fermi energy, Ep < (fc+/?A) 2 /2m < Ep+eV. The integral 
(|1.24p diverges logarithmically at k y — » 0. Assuming the same infrared cut- 
off in both subbands, k, we observe that the diverging In fc-contributions in 
the two subbands cancel each other, yielding in the linear (in V) response, 

f°° , eV /2A , v F + X\ , 

/ Szdx = —7- .„ In — (1.26) 

./_«, 2A(2tt) 2 \v f vf-XJ 17 

where Dp = \j2Ep~Jm is the Fermi velocity. Expanding this general result 
to the lowest non- vanishing order in X/vp we recover the net boundary 
polarization, Eq. (|1.5p . 

1.4.2. Evanescent modes 

We are now going to consider the edge spin density in the case of sharp 
(hard wall) edge potential. Since the sharp edge does not impose any 
finite length scale, the question arise, what would be the width of the edge 
spin distribution? Obvious candidate for that comes from the evanescent 
modes, 33 whose wave function do have an explicitly decaying component 
(~ exp( — m^/vpXx)) . 

The reflection at the sharp boundary mixes the two bulk subbands. 
Evanescent contributions in the upper subband appear when the reflecting 
states from the lower subband belong to the domain, m(vp — A) < k y < 
m(vp + A). Repeating the calculations leading to Eq. (|1.26|) but now for 
the evanescent domain only, we obtain, 

z 2A(2^) 2 \y vp l - 07^: J y ' 

Remarkably, the net evanescent contribution turns out to be much larger 
than the full result Eq. (|1.26p . This means that this contribution is largely 
cancelled by the contribution from the normal domain k y < m(vp — A). 
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Similar cancellation of all smooth long wavelength contributions takes place 
for the local spin density, as we discuss in the next section. 
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Fig. 1.4. Dependence of the local spin polarization, in units of eVm/8n 2 , on the dis- 
tance to the boundary for different values of spin-orbit coupling constant. Solid (red) 
line: X/vp = 0.1, dotted (blue) line: X/vp = 0.2, solid (black) line utilizes the approxi- 
mate formula Ql.300 for X/vp = 0.2. 



1.4.3. Hard wall, results 

Exact lengthy explicit expression for s z (x) in the case of sharp boundary 
potential, Fig. 11.41 was given in Ref. [35]. Here we show only the spectral 
density 



M) 



dxs z {x) sin qx, 



(1.28) 



which is given by a simple piecewisc continuous algebraic function defined 
in four domains. The surprising feature of the spectral distribution is its 
vanishing, s z (q) = 0, in the whole long- wavelength domain, < q < 2k + . In 
particular, this shows the exact cancelation between normal and evanescent 
modes. For larger values of q we obtain to the leading order in A, 



0, q < 2m(v F - A), 
eVq J —1, 2m(vF — X) < q < 2mvF, 
2mvp < q < 2m{vp + A), 
2m(vF + A) < q. 



(1.29) 

The plot of the spectral distribution is illustrated in Fig. 11.51 Remarkably, 
the net spin polarization (given by it" 1 J dqs z (q)/q) comes from the large-q 
tail (oc q~ 3 ) in the spectral density s z (q). 
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The approximate spin density may be written in a simple form (H = 1), 
eV 



s z {x) 



■ cos (2mvFx) sin (mXx). 



(1.30) 



2tt 2 VfX 

It is remarkable that the spin-orbit coupling constant enters via the period 
of beating only. 

As is evident from Eq. (|1.29[) one should speak about the spin accumu- 
lation at the hard wall with certain caution. The spin density in this case 
comes from quickly oscillating functions and the notion of spin accumu- 
lation should be understood in the same mathematical sense as the finite 
value of the integral 



dx sin x — lim 

J7— i-O 



dx 



' sin a; = 1. 



(1.31) 
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Fig. 1.5. Spectral distribution 111.291 1 of spin density in units of eV/invp for different 
values of spin-orbit coupling constant, X/vp = 0.1, and X/vp = 0.2. 



1.5. Smooth edge: Kinetic equation 

The second ballistic spin-Hall problem that allows analytical solution in- 
volves a smooth boundary potential, the exact condition to be presented 
below. Away from classical turning points (see the next section) spin ac- 
cumulation can be obtained with the help of the matrix kinetic equation. 
In addition this method allows to find a complete electron distribution 
along the way. In the semiclassical approximation the 2x2 matrix f p (x) 
generalizes the usual distribution function. The equation for the matrix 
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distribution function in the presence of a smooth confining potential U(x) 
and in the absence of disorder takes the form, 45 



Px f " 1 * d f 

~d x fp - A iySy, d X fpj + 2i\ [pySx ~ PxSy, fp] = ^ 



To the zeroth order in the potential gradient d x U the solution is trivial and 
reduces to the equilibrium Fermi-Dirac distributions for the two spin-split 
subbands (f3 — ±1), fp(x) = n F(^(p + /?A) 2 + U(x) — /it), with the local 
value of the potential U(x) determining the elevation of the bottom of the 
subbands. To this order the electron spin remains adiabatically within the 
plane of 2DEG during the course of electron motion in the potential U(x). 
The out-of-plane spin component arises from the non-adiabatic corrections 
that are first order in d x ll. The solution of kinetic equation (|1.32[) is rather 
straightforward and yields, 

/ p = \ £[1 + 20(n y s x - n x s y )]U + ^9*UJ2 (1-33) 

where n is the direction of the electron momentum. The local value of 
spin-Hall density s z — Tr s z J2 P fp is obtained by integrating Eq. (|1.33|) 
over excess electron states that originate in the left lead, namely over those 
with v y > and the energy within the interval [Ep, Ep + eV], 

'.(*) - I , ^ (---) [n F {E-eV)-n F {E)\. 

2W dx 2(2tt) 2 d\ A J y/2m[E-U{x)] \P+ P- J 

(1.34) 

Here p± are the momenta corresponding 

to a given energy: p±(x) = yj2m[E — U(x)] =F mX when both subbands 
are occupied, y/2m[E — U (x)] < mA; and p± — m\ qp y/2m[E — U(x)] 
when the upper subband is empty, y/2m[E — U(x)] > mA. Calculation of 
this integral in the linear order in eV yields 

; fj) - X " eV dU (135) 

z{ ' - {2^fm{v 2 F - A 2 )4 dx ■ ( ' 

The expression (| 1 . 35[) has a number of interesting features. The local 
spin polarization is proportional to the force exerted on electrons by the 
boundary. As long as the Fermi energy is well above the bottom of the 
bands, vp 3> A, spin accumulation is small and only second order in the 
spin-orbit coupling constant. When the bottom of the band is elevated 
high enough, vf < A, the local spin polarization increases dramatically. 
The vicinities of the two singularities at vp(x) — (classical turning point) 
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and v F (x) = A (degeneracy point) have to be studied by means beyond 
kinetic equation (jl.32p . 

In addition to a smooth classical spin distribution p.35[) quantum wig- 
gles in s z {x) are present whose magnitude is not necessarily small compared 
to s z (x). A convenient quantity (especially for numerical calculations) that 
averages out these wiggles is the "integrated spin density" 



Integrating Eq. (|1.35p over x yields smooth part of the integrated spin 
density 



This formula presents a generalization of Eq. (|1.26p for the case of a smooth 
variation of the confining potential. 

Note that the net spin polarization across the edge is again independent 
of the shape of the boundary potential, <S(oo) = — X 2 eV /12tt 2 v f , and is 
expressed via the bulk value of the Fermi velocity vf{oo). Here by x — oo 
we assume a point deep inside the 2DEG but yet far from its opposite edge. 
The latter has spin accumulation of the same absolute value and opposite 
sign. 

1.6. Smooth edge: singular spin dynamics near turning 



Scmiclassical Eqs. (|1.35ll.37|) do not offer the important information about 
the edge spin. First, these equations predict a singular spin density at 
vp(x) = and vp(x) = A. In addition, Eqs. (|1.35ll.37|) ignore any in- 
terference effects, which may be important for realistic boundary poten- 
tials. Both these problems may be addressed analytically, as it was done in 
Ref . [36] . In order to have more pedagogical discussion here we concentrate 
mostly on numerical results. Let us approximate the boundary by a linear 
potential 



with the constant force F. 

Figure [TTB1 shows the "raw" numerical data for s z (x) and S(x). The 
smoothness of the boundary implies that F <C m 2 X 3 /H. We see from 




(1.36) 




(1.37) 



points 



U(x) = -Fx, 



(1.38) 



December 23, 2009 15:2 World Scientific Review Volume - 9in x 6in Spin-Orbit-Book 



18 F. Author and S. Author 



14 
12 
10 

8 

N 

^4 

2 

-2 

20 40 60 80 

X 

Fig. 1.6. a). Spin density for a force strength F = 0.01m 2 A 3 /h (U = —Fx). [The curve 
offset both vertically and horizontally.] All coordinates are measured in units of h/mX. 
A smooth component of the density is hardly visible because of oscillating contributions, 
b). The smooth component is recovered in the integrated spin S(x) = s z dx' . Blue, 
green and red curves show the integrated spin for the force strength F = am 2 A 3 /h with 
a = 0.25, 0.05, 0.01, respectively. In all these cases we see three regions with different 
spin behavior. First, the spin density is the largest in narrow outer (x ~ 0) strip along 
the edge. This spin is compensated (and overcompensated) by the wide strip of negative 
smooth spin density. Finally, in the third strip (at x > 50 for the red curve a = 0.01) the 
smooth component of the density changes sign to positive again. The width of all three 
strips and the amount of accumulated spin, which in each strip is much larger than in 
Eq. 11. 5> . increase [formally unlimited] with decreasing slope of the boundary potential. 
To obtain the values of the spin density and accumulated spin one need to multiply the 
numbers in the figure by eVm/8hn 2 and eV/8\n 2 respectively. 

Fig. II. 61 a. how the rapid quantum oscillations make it hard to observe the 
mean value s z (x) Eq. (|1.35|) even for Fh/m 2 X 3 = 0.01. The smooth com- 
ponent is recovered in the integrated spin on Fig. II. 61 b even for relatively 
steep boundary Fh/m' 2 X 3 = 0.25. 

Explicit comparison between numerics and analytical expression 
Eq. (|1.37[) is made on Fig. 11.71 where we plot the rescaleed numerical results 
for different values of the slope of the boundary potential. 

According to Eq. (|1.37p we find two regions of different smooth spin 
behavior. First, within the strip where < Vp(x) < X spin density is 
negative (which is seen in a downward slope of the integrated density S(x) 
in Fig. II. 7[) . Farther away, s z {x) changes sign for Vp(x) > A, where both 
s z (x) and S(x) decrease gradually with increasing x. 
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Fig. 1.7. Integrated spin density S(x) = J^^Szdx for the potential U(x) = 
—ctm 2 \ 3 x/h in units of eV/8\ir 2 . The curves for a = 8,4,2,1 X 10 — 3 are shown in 
yellow, green, blue and red respectively. The horizontal coordinate is scaled differently 
for different curves, as x is measured in units of 10 3 a X h/mX. Narrow black lines stand 
for the semiclassical prediction, Eq. 11.261 . The logarithimc behavior, ~ log ct, of the 
dip at at U(x) = —m\ 2 /2 (x = 500) is clearly seen. Inset magnifies the region near the 
edge of 2DEG (asRiO). 

The most interesting is the behavior of spin at the borders of these 
regions, Dp = and vp = A. At Vf{x) = A the accumulated spin S(x) 
Eq. (|1.37| diverges logarithmically. This singularity originates from the ac- 
cumulation of classical turning points taking place when the conical cross- 
ing point in the spectrum of the Hamiltonian (|1.6p . see Fig. I1.3C passes 
through the Fermi energy. This singularity is regularized as S ~ \ogF 1 
according to Fig. 11.71 

Yet more peculiar is the behavior of both s z (x) and S (x) at the edge of 
2DEG, near the point where Vp(x) = 0. The smooth part of the accumu- 
lated spin, Eq. (|1.37p . has an infinite jump here (from Eq. (|1.37p it follows 
that <S(+0) = oo, while obviously <S(— 0) = 0). Development of such jump 
with decreasing slope of the potential is seen in the inset in Fig. 11.71 The 
jump in S(x) corresponds to the formation of a narrow strip with extremely 
large values of spin s z > along the border. This behavior will now be 
analyzed in more detail. 

Classical dynamics of electrons with Rashba spin-orbit interaction is 
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described by the effective Hamilton Eq. (|1.7[) . with the boundary potential 
approximated by the linear function U — —Fx. The family of classical 
trajectories generated by this Hamiltonian, shown in Fig. 11.81 demonstrate 
a number of unusual features. 

As seen from Fig. 11.81 those electrons from the lower subband that have 
\p y \ < mX, pass three turning points in the course of their motion in the x 
direction, corresponding to three solutions of the equation dH c ff /dp x = 0. 
Two of these turning points (those with p = mX) correspond to simultane- 
ous vanishing of both velocity components, v{x) — 0, the behavior generi- 
cally impossible in a 2DEG with the parabolic spectrum, H — p 2 /2m. 
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Fig. 1.8. Bottom: Family of classical trajectories at E = Ep for different values of p y , 
and m = A = F. Trajectories for both lower (red) and upper (green) spin-orbit split 
subbands arc shown. An example of a trajectory contributing to the peak in spin density 
at a; Ri 0, Eq. 111.401 1. is shown for "p x <C mA" (see the text). Top: The electron density 
p = ^\ipi + V , 2 1 /'2 for given longitudinal momentum and energy. Here Vl,2 are two 
eigenfunctions with p y = 0.2mA, Ep = 0, F = 0.003m 2 A 3 /?!, p in arbitrary units, x in 
units of h/m\. Three classical turning points can be seen. The interference of incoming 
and reflected waves in the upper subband causes smooth oscillations to the right of the 
inner turning point I (x > m\ 2 /2F). At the other turning points, II and III, the two 
kinds of oscillations are seen. Slow oscillations are caused by the interference of the 
incoming wave and the wave reflected at the turning point. Fast oscillations (wavelength 
~ h/mX) are due to the interference of distant (in time) segments of the same trajectory. 
Inset: Kinetic energy (arbitrary units) T±(p x ) = (p ± mA) 2 /2m for fixed p y = 0.2mA. 
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Analytical treatment 36 of the expectation value of the z-component of 
electron spin in the vicinity of turning point vf{x) = yields (x <C m\ 2 / F) 



where x = x(2Fm/h 2 ) 1 ^ 3 . In the asymptotic region x 3> {2Fm/h 2 )^ 1 / 3 
one can average over the oscillations of the Airy function. This allows us 
to recover the singular behavior of the smooth spin density (|1.35[) : (s z ) ~ 
x -3/2_ rj-v^g integral in Eq. (| 1 . 39|) features a logarithmic singularity at x — 0. 
With the logarithmic accuracy the (properly regularized) height of the peak 
of spin density is 36 



Striking feature of this result is that this maximal value is virtually inde- 
pendent of the strength of spin-orbit coupling or the shape of the boundary 
potential (up to a weak logarithmic factor) . 

We thus conclude that the nonequilibrium spin-Hall spin accumulation 
near a smooth boundary of 2DEG ballistic conductor with spin-orbit inter- 
action develops a narrow peak at the edge, with the width ~ (h 2 /mF) 1 / 3 
and height given by Eq. (|1.40[) . It is followed by a slow non-monotonic 
decay, as shown on Fig. 11.71 This smooth tail of spin density persists to 
much larger distances, > m\ 2 /F. The amount of spin accumulated in 
the peak (found as a maximum of the function S(x) = f x s z dx) equals 
<S max ~ 0MeV(m 2 /hF) 1 ^ 3 > 0.04eV/A, where in the last inequality we 
utilize the fact that F < m 2 X 3 /H. We see that the spin accumulated at the 
edge described by a semiclassical boundary potential is inversely propor- 
tional to the strength of spin-orbit interaction and becomes progressively 
larger for smoother slopes. This prediction can be used for experimental 
observation of spin-Hall effect in realistic two-dimensional electron systems. 

1.7. Spin-Hall effect on surfaces of topological insulators 

The method developed in Sections 11.41 and 11.51 is fully applicable to other 
chiral systems mentioned in the Introduction, as long as the correspond- 
ing terms in the Hamiltonians are linear in electron momentum. This is 
a good approximation in graphene and an acceptable one in Bi2Te3. In 
graphene, however, chiral structure arises from pseudospin (in sublattice 



l 




(1.39) 



o 




(1.40) 
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space). There is nothing straightforward about experimental detection of 
pseudospin in graphene. In addition, as explained in the Introduction, the 
case of graphene is mapped onto "Rashba-like" chirality for one cone within 
the first Brillouin zone and "Dresselhaus-like" chirality for the other one. 
As a result, the net pseudospin accumulation vanishes when contributions 
from the two cones are added together, cf. Eq. (|1.5p . 

The situation changes dramatically for topological insulators. First, 
the chirality originates from true spin, accessible by standard experimental 
techniques. Second, there is an odd number of Fermi valleys per (geometric) 
surface. This is in a sharp contrast to graphene which has two Dirac cones. 
Note that despite the equal number of "species" with opposite chiralities 
residing on the opposite surface, one can probe them individually ; by virtue 
of their spatial separation. 

Since the Fermi energy in topological insulators lies within the bandgap 
for 3D bulk electronic states, electron transport occurs only across the 
crystal surface. While not yet implemented in practice to our knowledge, it 
makes possible local control (via gate electrodes) of spatial distribution of 
electron density. We envisage the following experimental setup, see Fig. 11.91 
Metallic contacts are attached to the surface of topological insulator and 
drive dc electric current along y-direction. Gates are positioned (without 
direct contact) some distance above the surface. Application of electric 
potentials to the gates modulates position of electronic bands. Here, sim- 
ilarly to the rest of the paper, we assume that the topological insulator 
is disorder- free (ballistic). At the present time it is unclear how well this 
assumption is satisfied in contemporary samples, but undoubtedly sample 
quality is only going to improve in the near future. 

Effective low-energy Hamiltonian for Bi2Te3 46 (neglecting cubic terms) 



has a single Dirac point whose elevation is determined by the gate potential 
U(x). There are two cases of interest here: 

i) Potential U(x) is not very strong, so that the Dirac point always lies below 
the Fermi-energy (which is also normal situation for ungated material). 

ii) Potential is strong enough to lift the Dirac point above the Fermi level 
in some region on the surface, thus forming a p-n junction. 

We now present our general method before analyzing these two cases 
separately. The method extends the approach of Sec. 11.41 for the Hamil- 
tonian (|1.4ip and we present it here in detailed form for the reader's con- 
venience. Utilizing the fact that in a ballistic system k y is the integral of 



H = -iv(a x d y - o- y d x ) + U(x), 



(1.41) 



December 23, 2009 15:2 



World Scientific Review Volume - 9in x 6in 



Spin-Orbit-Book 



Spin-Hall Effect in Chiral Electron Systems 




Fig. 1.9. Proposed experimental setup for creating a n — p junction on the surface 
of topological insulator. Voltage V applied between the leads drives electric current in 
^-direction and creates a finite spin density along the junction (along the line x = 0). 



motion, we write the electron operators in the mixed representation, 



(1.42) 



(1.43) 



(1.44) 



Similarly, spin density in the mixed representation reads 

s(k y ,x) = (c{ y (x)ac ky (x)) 
and satisfies the equation 

dtSy (ky , ^c) = Ox Jx i^y > ^0 'ZvkyS z{ky , x) . 
Here J% is spin current operator (in general Jj = \ {piVj + VjCFi)), 

J«(k y ,x)=-v(cl y (x)c kv (x)). (1.45) 

In a steady state the left-hand side of Eq. (|1.43[) vanishes and we arrive at 

1 



2vk,, 



-d x J x (k y , x). 



(1.46) 
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This yields the net spin polarization across the profile of the gate potential 
U(x) 

s z {x)dx = -1^; 1^(^,00) - J*(ky,-oo)], (1-47) 



V 



expressing it via the values of spin current far away from the region where 
U{x) is varied, 

J%(k y , ±oo) = -vy]n(k x ,ky, ±oo), (1-48) 

with n{k x , k y , oo) denoting the corresponding distributions of electrons. 

The summation in Eqs. (jl.47jl - (jl.48p has to be performed over "uncom- 
pensated" states that originate from the left lead and belong to the energy 
interval Ep < E < Ep + eV, where V is the bias applied between elec- 
trical contacts. (The uncompensated states originate in the contact that 
has higher chemical potential and propagate towards the other contact.) 
The subsequent analysis will be performed separately for the two cases 
mentioned above. 

i) n-n junction. When gate potential U(x) is weaker than needed to 
elevate Dirac point above the Fermi-energy the electric current is carried by 
electrons only. Still, a smooth step-like potential U (x) create a "junction" 
between half planes with different values of the Fermi momenta, k^ and 
&r, respectively. 

Since dk x = dEk/k x v and dE — > eV the summation over k x is performed 
according to 

k m Ky^n k y 

(note that both k x > and k x < contribute to this expression). From 
Eq. (TOTD . 

^ dky k L dky k R 



r f \i ev ) f dk v k L f 

/ s z (x)dx = - — < / 

J-oo 27TU Jo k y k 2 _ k 2 Jo 



kl-kl JO h ^ k 2 R _ k 2 

(1.50) 

Regularizing these formally divergent integrals (by assuming the same lower 
cutoff which subsequently drops out from the difference of the two terms) 
we find the amount of net spin accumulated at the junction, 

s z {x)dx = ^\n(^f). (1.51) 
2nv V k R I 
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In addition, if the gate potential changes smoothly on the scale of 1/kp the 
local spin polarization can be written as, 

, , eV d x U(x) 

s ^ = ^dM- (L52) 

ii) p-n junction. In the real experiment the strength of the gate potential 
U(x) can be made significant enough to lift Dirac point above the Fermi 
energy over some region of the surface. For this setup we find 

s z (x)dx = 



DC 




(1.53) 

This result differs from Eq. (|1.50[) . and from the spin accumulation around 
vp(x) — A considered in section [L5l in that here the excess electrons to 
the left and to the right of the potential step have similar velocity, v y , 
but opposite momentum, p y , in the direction of current. In Eq. (|1.53p the 
excess electrons originate from the lead Ep + eV and thus have to have 
v y > 0. On the contrary, in the case of edge spin accumulation (Sec. 1 1 . 5j) 
the electrons arrive to the strip < vf(x) < A from the bulk 2DEG, as 
is shown in Fig. II. 8[ and thus must carry the bulk momentum p y > 0, 
but may have the "wrong" sign of velocity. As a result the spin density in 
Eq. p.53[) has the same sign on both sides of the n — p junction and two 
divergent logarithms add. Introducing a proper quantum mechanical cutoff 
to these integrals yields 



eV {4k p k n \ 2 1 dU 
s z {x)dx = - — In — — , g min = — — 
2nv V 9mm / hv dx 



(1.54) 



x=0 

Electrons with very small k y start to tunnel thorough the n — p junction, 
which would change our semiclassical predictions. The rude estimate of 
the quantum mechanical tunnelling exponent gives k = f " k x {x)dx = 
kyhv I (2dU I dx) . Requiring a small tunnelling probability, e~ K <C 1, gives 
the above cutoff q mui . 

We may now compare the strength of the Spin-Hall effect in semicon- 
ductor heterostructures Eq. (|1.5p with the result for topological insulators 
Eqs. (|1.51|1.54p . The main difference is that in conventional semiconductors 
the chiral electron states appear because of the weak spin-orbit interaction. 
Consequently the result Eq. (|1.5p acquires a small factor (A 2 — \ 2 D )/v F . 
There is now such suppression for the surface states in topological insula- 
torfl 

g Spin accumulation Eq. 11.51 and Eqs. 11. 5111. 541 1 depends also on the value of the 
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1.8. Summary 

In this paper we discussed several nonequilibrium spin-related phenomena 
occurring in 2-dimensional chiral electron systems. 

We were mostly interested in the effects which may be explained in 
terms of semiclassical electron motion in smooth potentials. In this case it is 
easy (at least theoretically) to produce the strongly in-plane spin polarized 
currents (Sec. II. 3| . 

The main part of the paper was devoted to calculation of the out of 
plane spin density (cr z ). It is the out of plane spin component, that is 
usually investigated in the experiment via the measuring of the optical 
Kerr rotation. 11 

Finally we discussed the mesoscopic spin-Hall effect in novel materials 
such as graphene and topological insulators. Experimental investigation 
of pseudospin-Hall effect might prove not to be easy. On the contrary, 
observation of spin-Hall a z accumulation at the surface of 3D topologi- 
cal insulator looks most appealing. While the net spin accumulation in 
semiconductor heterostructures is of the second order in small spin-orbit 
interaction, s ~ X 2 /vp (|1.5p . there is no such restriction for topological 
insulator, s ~ 1/v (jl.51ll.54[) . As a result, the spin accumulation in the 
latter case should exceed that possible in traditional 2-dimensional electron 
gas by orders of magnitude. 
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